#  destripe v 0.4 - Removes stripes from tiles of geophysical data.
#                   AKA Zero mean traverse.
#  Copyright (C) 2010 John Donovan - mersey.viking@gmail.com

#  Licensed under the Apache License, Version 2.0 (the "License");
#  you may not use this file except in compliance with the License.
#  You may obtain a copy of the License at

#  http://www.apache.org/licenses/LICENSE-2.0

#  Unless required by applicable law or agreed to in writing, software
#  distributed under the License is distributed on an "AS IS" BASIS,
#  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
#  See the License for the specific language governing permissions and
#  limitations under the License.

def destripe(image, horizontal=True):
    for tile in image.tiles():
        tileMeanValues = None

        if horizontal:
            tileMeanValues = tile.mean(axis=1)
            tileMeanValues.shape = (tile.shape[0], 1)
        else:
            tileMeanValues = tile.mean(axis=0)
            tileMeanValues.shape = (1, tile.shape[1])

        tile -= tileMeanValues

